Experimental Study of Noise-induced Phase Synchronization in Vertical-cavity Lasers 
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We report the experimental evidence of noise-induced phase synchronization in a vertical cavity 
laser. The polarized laser emission is entrained with the input periodic pump modulation when an 
optimal amount of white, gaussian noise is applied. We characterize the phenomenon, evaluating 
the average frequency of the output signal and the diffusion coefficient of the phase difference 
variable. Their values are roughly independent on different waveforms of periodic input, provided 
that a simple condition for the amplitudes is satisfied. The experimental results are compared with 
numerical simulations of a Langevin model. 

PACS numbers: 05.40.-a,42.55.Px 



The phenomenon of Phase Synchronization (PS) has 
been the subject of extensive investigation during the 
last years (see e.g. 0] for a recent review). In the 
purely deterministic case, two systems with instanta- 
neous phases 3>i(i) and $2(i) are locked if, for two in- 
tegers n, m, n$i — m<i>2 remains bounded for all times. 
In the presence of unbounded fluctuations, this condition 
on the phase difference cannot hold for all times. Suffi- 
ciently large noise will eventually cause phase slipagge 
and PS can thus only be effective. For example, for weak 
periodic forcing of an autonomous stochastic oscillator, 
the dynamics of the instantaneous phase difference <t> is 
ruled by the Adler equation 



(i) 



where £ is a Gaussian white noise. Eq. Q describes 
the motion of an overdamped Brownian particle in the 
tilted potential —Acf> — A s cos0. Therefore, if the fre- 
quency mismatch A is smaller than the synchronization 
bandwidth A s , locking can occur. However, from time to 
time fluctuations suffice to kick <j> out of the well yielding 
a phase slip. Recent observations of noise-enhanced PS 
have been reported in more complex situations, where 
chaos occours in absence of noise jjj Q. In this con- 
text, dynamical features may be significantly altered and 
a richer variety of phenomena arise Q . 

For noise-driven bistable systems, PS is strictly related 
to the phenomenon of Stochastic Resonance (SR), i.e., 
the enhancement in the response to a weak signal super- 
imposed to a stochastic input (see e.g. 7;]). The very 
fact that noise can induce synchronized jumps between 
the two stable output states has been recently interpreted 
as a noise-induced PS. A first observation of PS in an 
electronic circuit has been reported in || , analyzing the 
output mean frequency as a function of added noise. Nu- 
merical simulations of Langevin models 0, llfj| showed a 



similar behavior and allowed for a more refined descrip- 
tion of PS in terms of the effective phase diffusion co- 
efficient 6]. An analytic theory, performed in the case 
of a periodic and aperiodic square wave modulation, was 
presented by Freund et al. in where the explicit ex- 
pressions for the mean output frequency and the effective 
diffusion coefficient were derived from a Master Equation 
for the stochastic variable cj). 

In this work, we present a detailed experimental in- 
vestigation of noise-induced PS regimes in a a bistable 
optical system. Our study is motivated by the necessity 
to perform a throughful test of the statistical indicators 
of PS in reliable, well controllable laboratory conditions. 
The effect of both sinusoidal and square-wave input pe- 
riodic modulations are considered; the latter case allows 
to discuss a general framework for a comparison with the 
analytic results of Ref . . 

Our experimental system is a pump-modulated Verti- 
cal Cavity Surface Emitting Laser (VCSEL). The VC- 
SEL is a semiconductor laser, with a symmetrical cav- 
ity allowing for the emission on two linear, perpendicu- 
lar polarizations selected by the crystal axis directions 
(see e.g. [Tij ;. The emission symmetry, usually bro- 
ken by impurities or optical inhomogeneties leading to a 
single-polarization emission, may be restored for partic- 
ular choices of the pump current. In such a case, the sys- 
tem is bistable and exhibits noise-driven, random jumps 
between the two polarizations ruled by the Kramers' law 
[l3| . The experimental evidence of SR 01 an d of binary 
aperiodic SR 0] in a VCSEL has been recently reported. 

In our setup, we employ a VCSEL lasing at 850 nm, 
thermally stabilized (better than 1 mK) and with a 
carefully controlled pump current. The overall stabil- 
ity allows for long-time measurements, even in presence 
of critical behaviors. A linear polarization direction in 
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FIG. 1: The Spectral Power Amplification versus applied 
noise amplitude for different amplitudes of SIN modulation, 
Asin = 50, 100, 150, 200, 250 mV (top to bottom). The inset 
is a magnification of the resonance region. 



the laser emission is selected using a polarizer and an 
half-wave plate. The laser intensity is monitored by an 
avalanche detector and the signal is recorded by a digital 
scope. An optical isolator prevents from optical feedback 
effects to occur. The signals from a 10 MHz-bandwidth 
white-noise generator and a sinusoidal (SIN) or a square- 
wave (SQR) oscillator are summed and coupled into the 
laser by means of a bias-tee. The period of the input sig- 
nals was chosen to be T = 5 /j,s (200 kHz) for both wave- 
forms. The measurements are performed varying the in- 
tensity D of the applied noise as well as the amplitude 
A of the input modulation. For all the measurements 
reported henceforth, the A values are subthreshold, i.e. 
small enough to avoid output transitions every modula- 
tion period (in absence of added noise). The amplitude 
Asqr is set to be equal to the root mean square value 
of Asin, i-e. Asqr — Asin/V%- For every choice of the 
parameters, a sequence of 16 million points is acquired, 
corresponding to 500 samples/period for a total of 32,000 
periods. Such large data sets allow to reduce statistical 
errors to a level of accuracy comparable with current nu- 
merical simulations of theoretical models. 

To illustrate the phenomenon of SR, we first report 
in Fig. the response of the system for different values 
of the amplitude of the SIN modulation, as evaluated 
by means of the spectral power amplification SPA = 
4 |X OMt (fi)| 2 /A| JAr 0, where X OMt is the Fourier spec- 
trum of the output signal evaluated at the driving fre- 
quency f2 = 2ir/T. A well defined peak of SPA for an 
optimal value of the applied noise is observed, yielding 
the typical signature of SR • Increasing the input am- 
plitude, a decrease in the response as well as a shift to- 
wards lower noise values of the peak occurs, as already 
reported from numerical simulations and observations in 
nonlinear circuits 1171 . 




FIG. 2: The phase difference cj> for SIN (left) and SQR (right) 
inputs of amplitude Asin = 250 mV . Noise levels are 0.0121, 
0.0169, 0.0225 Vrms (bottom to top). The insets are the time 
series for low-amplitude (Asin ~ 100 mV) and noise level 
close to resonance (0.0196 V^ Mg ). 

The peak in the SPA indicates the improvement of the 
response of the system at the modulation frequency, thus 
suggesting that a PS regime is achieved. To investigate 
such a possibility we define the input-output phase dif- 
ference </>(t) = $out(*) ^ $in(*) where $ in (<) = Ot. The 
instantaneous phase of a real signal x(t) can be defined 
in a general way employing the concept of the analytic 
signal [lj . This is accomplished by considering the com- 
plex signal z(t) = x(t) + iy(t), where y = fix is the usual 
Hilbert transform of x(t). The phase is than given by the 
argument of z. For bistable systems an alternative and 
easier definition can be given based on the knowledge of 
the sequence of switching times t n only Q]; $ out can be 
defined either as a piecewise-constant function increas- 
ing by steps of ir at each jump time t n of the output or 
by linear interpolation. In the analysis of experimental 
data we checked that different definitions yield the same 
results, up to statistical accuracy. 

In Fig. [5] we report the temporal evolution of 0, as ob- 
tained from the Hilbert transform method, for both the 
SIN and SQR waveforms and different values of the added 
noise. In both cases, a clear signature of PS is found at a 
well defined (resonant) noise power corresponding to the 
peak in the SPA. As shown in the Figure, for a noise be- 
low (resp. above) the resonant value the phase difference 
drifts towards lower (resp. higher) values. The phase dif- 
ference at resonance mantains an almost constant value, 
rarely interrupted by jumps of integer multiples of 2tt 
(phase slips). For comparison, in the insets of Figs. it 
is shown the temporal behavior of 4> at resonance for a 
smaller value of the input signal amplitudes; the phase 
slips are here more frequent, yielding a faster diffusive 
motion. 

To quantify the observed behaviors, it is useful to in- 
troduce the average frequency (w ut) — ($out) and the 
effective diffusion coefficient V c g = \4? [((j) 2 ) — (<j>) 2 ] 
0. The measured values of (w ou t) and T> e g are reported 
in Fig. |3 as a function of the noise intensity for different 
values of the input signal amplitude A. Upon increas- 
ing A the region of locking, i.e. the range of noise val- 
ues for which (ui ou t) — widens while 2? c ff develops a 
pronouced dip at the resonant noise value. It is worth 
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FIG. 3: Average output frequency (a) and diffusion coef- 
ficient (b) for sinusoidal and square-wave inputs (full and 
open symbols respectively). Input amplitudes are As in = 
50, 150, 250m V (circles, squares, diamonds, respectively). 



noticing that, while the output frequency equals the in- 
put frequency at a well defined value D = D* , increasing 
A the abscissa of the minimum of T> c g shifts towards 
lower noise amplitudes. Moreover, such abscissa corre- 
sponds, within the experimental accuracy, to the abscissa 
of the maximum of the SPA. Such result is in agreement 
with the analytic predictions of Ref. 11]. 

A related way to illustrate the enhacement of phase co- 
herence is to plot the average duration of locking episodes 
{Tiock) (i-e. the average time between consecutive phase 
slips) as a function of the noise 0] . As seen in Fig. 01 
a marked maximum signals strong input-output correla- 
tions, namely a decreasing rate of phase slips (see again 
Fig. 0). The minimum value of the diffusion constant 
2? e ff and the maximal value of (Ti oc k) & re indicators of 
the quality of PS. It is therefore important to estimate 
their dependence on the experimental parameters. On 
the other hand, the two are not independent: as argued 
in Ref. 0] i the very definition of T> c g yields 
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D*, we expect (cb) ~ 0, and therefore 7r 2 /22? ff ~ 
fe). A shown in the inset of Fig. 0] this estimate is 
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FIG. 4: Average duration of locking episodes for SIN and 
SQR inputs (circles and squares respectively), As in = 
50, 100, 150,200, 250 mV (bottom to top). The inset shows 



(Tiock)/T and ir 2 /2V cB T for D 
A/D (see text). 



D* = 0.0196 Vr MS ) versus 



waveforms. 

A further relevant information is that (Ti oc k) increases 
exponentially with A/D, confirming that PS rapidly be- 
comes effective upon increasing the driving amplitude. 
A simple argument to understand this scaling is the fol- 
lowing. Consider the case of a two-level system with a 
slowly modulated barrier AV(t) = AV ± B, where the 
sign changes every semiperiod T/2. The two Kramers' 
rates ax^ — tk exp(^fB/D) give the probability per unit 
time to jump from one level to the other. During the 
locking period, the system stays synchronized with the 
modulation, likely jumping (almost) every semiperiod. 
But the probability per unit time for the system to jump 
less then (or more than) one time per semiperiod (re- 
sulting in a phase slip) is approximatively a\ for a sufli- 
cently large B/D. This accounts for the observed scal- 
ing if B is proportional to A. The same conclusion can 
be drawn from Eq. 0J. Close to resonance (A ~ 0) 
2?eff oc exp(— const. A s /D) 6] and the exponential in- 
crease then follows from the fact that A s oc A for small 
amplitudes. 

The experimental data are in qualitative ag reement 
with both the numerical 0] and analytical |llj studies 
of noise-induced PS. This originates from the fact that, 
in the SR regime, the polarization dynamics of our VC- 
SEL can be succesfully modeled by the phenomenological 
Langevin equation 

x = -V'(x)+Af(t)+r) , (3) 

where V is a double-well potential, / is the modulation 
waveform (of frequency Q) and 77 is a white, Gaussian 
noise such that (n(t)n(t')) = 2D5(t - t'). Here the di- 
mensionless variable x represents the polarized laser in- 
tensity. To emphasize the correspondence with the ex- 
periment we use here the same symbols A, D, Q as above. 
However, one should be aware that this correspondence 
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cannot be established a priori since the model is only 
an effective description. A direct quantitative compari- 
son is feasible only by means of a careful calibration of 
the physical parameters of model |yj on the experimen- 
tal data. This task has been accomplished in slightly 
different experimental conditions Q, and goes beyond 
the scope of the present work. Here, we limit ourselves to 
illustrate the results of simulations for the model poten- 
tial V(x) — x 2 (x 2 — 2) using a second-order stochastic 
Runge-Kutta algorithm. The outcomes for moderate val- 
ues of A (A < AV) shown in Fig. [S] compare well with the 
experimental ones (see again Figs.0). Furthermore, the 
exponential scaling T> c ff oc exp(— const. A/D) is found, in 
agreement with experimental findings. 

Let us now address the issue of how PS is affected 
by different choices of the driving. Figs. indeed show 
that both (wout) and T> c g are approximatively indepen- 
dent on the modulation waveform, provided that A$qr = 
As i n I a/2. The only relevant parameter is, at least in this 
range of amplitudes, the modulation frequency. To check 
the generality of this observation we performed simula- 
tion of Eq. @ with three different waveforms /, namely, 
a sinusoidal, a square-wave and a sawtooth (SAW) with 
the same RMS value, Asqr = Asin/V^ = Asaw/V^- 
As shown in Figs. [SJ the curves are almost overlapped in- 
dicating that for moderate amplitudes the quality of PS 
is largely independent of the details of the driving. To 
understand this behaviour, we compared the Kramers' 
rates for the different waveforms averaged on the modu- 
lation semiperiod. Upon changing A/D in the range of 
interest, we indeed found that they are almost indepen- 
dent on /. Therefore, in the adiabatic limit a comparison 
with the theory can be worked out and is currently 
in progress. 

In conclusion, we reported a detailed experimental ev- 
idence of noise-induced PS in a VCSEL. Phase entrain- 
ment is achieved by different input periodic signals, with 
the same value for the measured statistical indicators as 
long as the RMS of the signal is the same. The numerical 
analysis of a Langevin model reproduce the features ob- 
served experimentally, included the scaling of the locking 
times at resonance with the amplitude of the modulation. 
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